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Summary. Observations suggest that many, if not all, post AGB systems evolve 
through an aspherical outflow phase. Such outflows require a sufficient engine rota- 
tional energy which binaries can provide. Via common envelope evolution, binaries 
can directly eject equatorial outflows or produce poloidal outflows from magnetized 
accretion dislcs around the primary or secondary. We discuss how accretion driven 
magnetohydrodynamic outflow models all make similar predictions for the outflow 
^ . power and speed and we distinguish between the launch vs. propagation regimes 

^ ' of such outflows. We suggest that the high velocity bipolar outflows observed in 

C/3 , planetary nebulae (PNe) and the lower velocity but higher power bipolar outflows 

^ , ■ observed in pre-PNe (pPNe) are kinematically consistent with time dependent ac- 

cretion onto a white dwarf (WD) within a depleting envelope. Since the WD primary 
, core is always present in all post- AGB systems, accretion onto this core is potentially 

^ ' common. Previous work has focused on core accretion from sub-stellar companions, 

0^ ' but low mass stellar companions may be more important, and further work is needed. 

a^ ■ 

Key words: stars: AGB and post-AGB; (stars;) binaries: general; accretion, accre- 
tion disks; magnetic fields; stars: winds, outflows 



o 

O 



1 Introduction: Kinematics of pPNe and PNe 



Understanding the origin of asymmetric outflows in PNe and pPNe requires 
. feedback between observers, specific object modelers, and paradigm-seeking, 

^ order of magnitude theorists. Here we behave as the latter. 

Generally, pPNe exhibit a combination of a fast bipolar outflow embedded 
within a slow spherically symmetric wind from the AGB star [1]. Presently, 
the data do not rule out all pPNe having gone through a strongly asymmetric 
outflow and all PNe having gone through an asymmetric pNE phase. The 
symmetry of PNe would then correlate with age and the evolution from the 
pPNe to PNe could reflect a time evolution of the same physical mechanism 
that produces asymmetry on small scales but leads to a a nearly spherical 
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structure on large scale as supersonic motions damp. While AGB stars produce 

spherically symmetric outflows, pPNc asymmetry arises within < 100 yr [1, 2]. 

For pPNe [1], each fast wind has a typical age At ~ 10^ — lO^yr, speed 
~ 50km/s, mass Mf ~ 0.5Mq, outflow rate, Mf ~ 5 x lO~^M0/yr, momen- 
tum 77 ^ 5 X lO'^^g.cm/s, and mechanical luminosity L^j > 8 x 10''"''crg/s (can 
be as high as lO^^erg/s). The slow pPNe wind has an age At ^ 6 x lO'^yr, a 
speed Vw ~ 20km/s a mass Ms ~ 0.5Mq, outflow rate, Ms ~ lO~*M0/yr, mo- 
mentum Us ^ 2 X lO'^^g cm/s, and mechanical luminosity Lm.s ~ 10^'*crg/s. 

For PNe observations suggest e.g. [3] an age At ~ 10'*yr a slow wind of 
speed Vs ~ 30km/s of mass Ms ~ O.IM©, outflow rate, Ms ~ 10~^M©/yr, 
momentum Ug ^ 6 x lO^^g cm/s, and mechanical luminosity Lm.s ~ 3 x 
lO^^erg/s. PNe have fast winds of speed t;/ ~ 2000km/s, mass M / ~ 10~^Mq, 
outflow rate, Mf ~ lO~*M0/yr, momentum 77/ ~ 4 x 10'''g cm/s, and me- 
chanical luminosity L^.f ^ 1-3 x lO'^^crg/s. 

The pPNe phase demands the most power. The linear momenta of fast 
bipolar pPNe outflows seems too large for radiation driving [1] . This motivates 
the need to tap rotational energy which binaries can provide. The high fraction 
of binaries attributed to p/PNe has led to speculation that all asymmetric PNe 
and pPNe involve binaries [4, 5] . 

How is the rotational energy and angular momentum converted into out- 
flows? Differential rotation supplied by binaries can amplify magnetic fields 
which can in turn produce accretion powered bipolar jets. Accretion powered, 
magnetically mediated jets, arc scicmingly iibiquitous in astrophysics and can 
accommodate the high momentum pPNe demands [6]. 

2 Which Binary Accretion Scenario? 

In a common envelope (CE) [7] the companion drags on the envelope of the 
primary, transferring angular momentum and kinetic energy. Provided that 
the envelope cooling time scale exceeds the energy supply time from the in- 
spiraling companion, a converted fraction a > 0.1 of its loss in gravitational 
energy can spin up the envelope and unbind it. For a fixed envelope mass, 
a lower mass companion must fall deeper to unbind the envelope. Fig. 1 [8] 
shows aspects of CE for low mass companions in an AGB envelope. When the 
Roche radius is reached by the in-spiraling companion, accretion can occur. 

Ref [9] discusses the secondaries for which accretion disks will form around 
the primary core. Brown dwarf (BD) (0.003Mq < Merit ~ O.OTM©) radu 
increase with decreasing mass unlike their Roche radii which decrease with 
decreasing mass. Such objects will unstably lose mass. Since the circulariza- 
tion radius lies outside of the primary's core, a disk can form within a few 
orbit times. This contrasts the M > Merit case for which the stellar radius 
decreases with decreasing mass and more strongly so than the Roche radius. 
Supercritical companions present a circularization radius within the primary's 
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Fig. 1. From [8]: For a 3Mq young AGB star (left) and inter-pulse AGB star 
(right), the solid lines show the radius where a companion must in-spiral to unbind 
the envelope for different values of the drag parameter a. The dash-dotted line shows 
the core-envelope boundary. The long-dashed line marks the tidal shredding radius 
and the short dashed line is where the companion first fills its Roche lobe. 



core. Material leaving the secondary would then initially spiral swiftly into the 
primary rather than orbit quasi-stably. 

Ref.[9] focuses on BDs (O.OO3M3 < M2 < Merit), but accretion ft'om 
low mass stars and planets warrant further study: Fig.l (left) shows that for 
low values of the (uncertain) drag parameter a, some low mass stars could 
tidally shred into a disk before reaching the radius where they unbind the 
envelope and subsequent in-spiraling is potentially slowed. Fig. 1 also shows 
that planets will shred into a disk upon inspiral and for large a in the inter- 
pulsc AGB, even large planets might unbind the envelope. Gcnciralizations 
must incorporate the fact that the companion size can be of order the tidal 
shredding or Roche lobe radius, and the structural evolution of the primary 
as the secondary in-spirals. The latter is important for radii both exterior and 
interior to the companion: How fast do further inspiral and angular momentum 
transfer subsequently occur once the outer envelope is ejected? 

If a stellar companion unbinds the envelope at a radius with Roche lobe 
overflow but without tidal shredding, sustained accretion requires the binary 
to lose angular momentum. As the secondary fills its Roche lobe, it will drag on 
any residual inner envelope material. This could transfer the needed angular 
momentum, but too much drag could prevent the companion from forming a 
quasi-stable Keplerian disk. Material that initially leaves the companion for 
the primary also carries magnetic fields and thus will magnetically link the 
primary and secondary. Even though the circularization radius for M2 > Merit 
is inside the core, the in-spiraling material still incurs differential rotation. 
Magnetic fields can then be amplified and can also release angular momentum. 

Much of the initial accreted energy in the M2 > Merit case would be 
released upon material impact to stellar surface. This could produce dwarf 
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novae type bursts shrouded by the stellar envelope. However, if the right 

range of mass and angular momentum are transferred in this initial accretion 
phase to (1) drop M2 below Merit, (2) keep M2 filling its Roche lobe, and 
(3) leave enough angular momentum to form a Keplerian disk, then accretion 
could proceed as for the initial M2 < Merit case of Ref.[9]. Rcf. [10] suggests 
a ~ 0(100) year delay in the time scale between ejection of cicumbinary dust 
tori and presence of jets in pPNe. If the companion initially has M2 > Merit, 
the fraction of the envelope ejected as it spirals in might become the dust torus, 
and the delay before the jet could be the time it takes for the companion to 
lose enough mass to move the circularization radius outside the core. 

The role of low mass stellar and planetary companions is important be- 
cause although ~ 16% of nearby (< 50pc) sun-like (late F to early K) stars 
have companions more massive than Jupiter at < 3AU, 11% are stars, 5% are 
planets, and 0% arc BDs [11]. If this ^ 16% were crudely taken as the fraction 
of binary induced, and thus asymmetric p/PNe, this is within a factor of 2 of 
estimates of the total fraction of low mass stars incurring any PNe [5]. 

It is also possible for accretion disks to form around the secondary [12, 13, 
14]. The accretion rate inferred from the Bondi wind accretion formulae is 

M, \mJ [1 + (i,/x;,)2]3/2' ^ ) 

where v is the orbital speed of the secondary and Vw is the slow wind speed 
from the primary. In general, for M2 < Mi, reasonable parameters provide an 
accretion rate compatible with the luminosities rcriuircd at the PNe stage if 
the companion is a main scqTience star. However, the outflow velocities of the 
fast wind in PNe would require a WD companion, as seen in the next section. 
Complcmentarily, the luminosities of the fast wind pPNe outflows would also 
require a WD companion. 



3 Accretion Disk Outflows in pPNe and PNe: 

Accretion disk outflows have a mechanical luminosity of order [6] 

=4.5 X lO^^e-i— , (2) 

where e is the dimensionless efficiency of conversion from accretion to outflow, 
Ri is the inner disk radius. Ma is the accretion rate, G is Newton's constant, 
and M* is the central stellar mass. In the last term 10 = iJi/lO^^cm, e_i = 
e/0.1 and Ma,-4 = Ma/lO-'^Mg/yr. For an MHD outflow, Eq. (2) equals the 
Poynting flux at the launch surface. 

For propagation into a region of negligible inertia, the asymptotic outflow 
speed is ~ f2rA [15] where Q is the angular speed of field anchor point and 



Pre-Planetary/Planetary Nebulae Engine Principles and Paradigms 



5 



rA is the radius where the poloidal outflow speed equals the Alfven speed. 
This product is typically of order 1-3 times the escape speed of the inner most 
radius of the disk and is thus at least 

Vout ~ Vesc = 1600 — ^ km/s. (3) 

Time dependent accretion outflows described with Eqs. (2) and (3) are 
consistent with the high pPNe outflow mechanical luminosity and the fast 
PNe wind speed of Sec. 1 when M*/i?* corresponds to a WD (for accretion 
onto primary M* = Mi, or secondary M* = M2) The outflow speed from 
(3) does not depend explicitly on Ma (only implicitly and weakly via Ri 
which can exceed the stellar radius for a strong field) and is thus largely time 
independent. However, this speed depends on the inertia of material blocking 
the outflow. Conservation of momentum gives 

^obs = f ~ 80km/s, (4) 

where M//Mq = 3.3 x lO-^e-iMao.-s r-^/^'dT is the mass in one of 
the fast collimated outflows, Me„„ is the envelope mass, fo ~ 0.2 is the solid 
angle fraction intercepted by the collimated outflow, and r = t/lyr is used to 
incorporate Ma oc t~^^^ of Ref. [9]. The numbers have been scaled to the pPNe 
case so that here M^nv » Mf and for an envelope of mass 2Mq, giving the 
intercepted mass of 0.2 Mq for fa = 0.2. Eq. (4) then implies Vobs = 40km/s, 
within the range observed [3]. 

Eq. (4) represents the observed speed of the fast when blocked and loaded 
by the envelope. By the end of the pPNc phase, the envelope is quite ex- 
tended, reducing the optical depth and revealing material moving at the "free 
streaming" fast wind speed. Assuming a dust-to-gass mass ratio of 1/100 and 
micron sized grains of density of 2g/cm^, the optical depth from dust is 

^ ^2J71FTW^) (10^) {1^^ ' 



scaled for PNe. For pPNc, the density increases by a factor > 10'' and R 
is down by a factor of 10, so Td > 2.5, optically thick. The different optical 
depths of pPNe and PNe can thus explain why observed PNe fast winds can 
have Vf > 1600km/s, whilst those of pPNc have vj < lOOkm/s. 

Keeping in mind our discussion of low mass stellar companions in Sec 2, 
note that Ref. [9] considers a companion of mass M2 ~ O.O3M0 < Merit and 
a Shakura-Simyacv viscosity parameter ass ~ 0.01, for which the accretion 
rate then decays as Ma ~ 1-6 x lO~^t~^/''M0/yr.. Using this in (2) with 
e = 0.1 for t = 100 yr with Ri = 2 x lO^cm and Mi = 0.6Mq gives Lmj ~ 
4.3 X 10"^^(t/lyr)~'''/^. This provides the needed power demands of Sec. 1 for 
pPNe after 1000 yr and for PNe after lO'' yr. 
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Because the surface density evolves, the gas opacity evolves from Thom- 
son to Kramer's after ^ 100 yr [9], and the height to radius ratio decreases 
substantially as the disk cools. A fixed a^g ~ 0.01 is self-consistent with the 
time evolving accretion and power above. 

4 Launch versus Propagation Regions for MHD Outflows 

We refer to the "launch" region [16] of MHD outflows as the region where 
the magnetic force and energy dominates the flow and thermal energy. This 
extends to a height typically no greater than Zc ~ 10 — 50i?i, where Ri is 
the innermost radial scale of the engine, (e.g. the inner radius of an accretion 
disk). In the launch region the bulk flow is accelerating but is sub-Alfenic until 
reaching Zc- The "propagation" region describes z > Zc where the poloidal flow 
speed exceeds the Alfen speed, eventually approaching its asymptotic speed. 
Presently, only the propagation region is observationally spatially resolved. 
With this distinction, we describe 3 classes of MHD outflow related work. 

Launch to Propagation: e.g. [17, 18, 19]; Here the magnetic fleld is imposed 
to be ordered on a scale at least as large as the anchoring rotator. The base 
of the rotator is typically a boundary condition. Calculations address how 
material is accelerated by the combination of centrifugal flinging of material 
along quasi-rigid field lines and/or a vertical gradient in the magnetic pressure. 
Both a poloidal and toroidal fleld component at the base are required. The 
flow can be coUimatcd and supersonic by hoop stresses before upon reaching 
the propagation region provided that there is an ambient pressure to coUimate 
the magnetic field. Simulations can cover from the base of the launch region 
to a scale typically < 100 times the engine scale. 

Asymptotic Propagation: e.g. [20, 21]; Here a collimated jet is injected on 
small scales, and the subsequent propagation and shaping by the ambient 
medium of specified magnetic and thermal properties is studied. The jet and 
ambient medium parameters are varied to assess what conditions can produce 
the observed asymptotic, morphologies. These simulations do not address the 
magnetic field origin or acceleration mechanism. 

Field Origin to Launch: e.g. [6, 22]; The previous two categories do not 
address where the dynamically important large scale fields comes from in the 
first place. Accretion of flux may be difficult in a turbulent disk, but these large 
scale flelds can plausibly be produced by a combination of a flow dominated 
helical dynamo inside of the rotator, followed by magnetic buoyancy, and a 
magnetically dominated helical dynamo relaxation in the corona [23]. The 
latter opens up structures to the large scales needed to drive jets, much as 
solar coronal loops open to create solar coronal holes. This category of work 
focuses on the field origin, with kinematic estimates of the subsequent launch, 
but not the dynamical launch itself. 

It is a frontier to couple the above categories, and each has limitations 
when separated from the other two. For example, a self-consistently grown 
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strong field dynamically mediating the launch need not necessarily imply a 
strong magnetic influence in the propagation region: If the launch region pro- 
duces a super-magnetosonic coUimated outflow, then even subsequent ballistic 
propagation into the propagation region would still emerge as collimated. In 
addition, a super-magnetosonic outflow could become turbulent and the tur- 
bulence can amplify small scale magnetic energy to > 10% of equipartition 
with this turbulence. This can be a substantial fraction of the initial bulk flow 
energy and such a flcld would be responding to the flow, not the reverse. 

To summarize: the physics of the launch region (< lAU) involves: (1) Ori- 
gin of large scale magnetic fields, field buoyancy to coronae, field relaxation 
into larger coronal structures, (2) physics of centrifugal + magnetic accelera- 
tion of material from small to super- Alfvenic speeds, or Poynting flux driven 
bursts of acceleration, (3) criteria for steady or bursty jets, and (4) assessment 
of the extent of Poynting flux domination. 

The physics of the propagation region (>> lAU observationally resolved) 
involves such issues as: (1) Propagation, instability formation, and sustenance 
of coUimation in as a function of internal vs. external density and strength 
of magnetic flelds, (2) bow shocks, cocoon physics, particle acceleration, (3) 
effect of cooling on morphology, and (4) interaction with ambient media 

5 Toward Connections between Theory and Observation 

Spatially resolving the launch region and measuring the magnetic field strength 
and geometry therein would be the gold standard for directly evaluating the 
role of magnetic fields in producing asymmetric; p/PNe. Measurements of fields 
in the propagation region provide primarily indirect evidence, though the de- 
tection of relatively strong fields there is particularly significant. [2]. 

Whether binaries supply needed rotational energy to amplify jet-mediating 
fields in accretion disks is a fundamental question. The basic wind kinematics 
are roughly consistent with accretion onto a WD, suggesting the importance 
of accretion onto the primary. 

Coupling accretion disk physics to large scale magnetic field production, 
to jet launch, and jet propagation in a unified theory is work in progress that 
spans several subfields of theoretical astrophysics, let alone the specific ap- 
plication to p/PNe. However some predictions/trends can be studied: (1) It 
should be possible to evaluate the kinematic constraints/predictions of Sec 3 
in more detail and compare the distribution of inferred fast outflow speeds 
to what would be expected from known binary statistics of low mass main 
sequence stars. This would serve to help determine the commonality of accre- 
tion onto the companion vs. the primary. (2) CE models would predict mostly 
Oxygen rich rather than Carbon rich post AGB systems because typically, the 
companion ejects the envelope on time scales of order years, right from the 
beginning of the AGB phase when the envelope expands, so the AGB star 
would not have had a chance to reach the Carbon dredge up. (3) Crystalline 
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dust in post-AGB systems can be produced if a binary induced spiral shock 
anneals silicates [24]. Is this universal? (4) CE evolution would predict equa- 
torial outflows from companion inspiral that precedes any accretion driven 
poloidal jet. A delay is observed [10] but more work is needed to predict the 
delay time scale. (5) Arc the geometry and composition of dust tori around 
post-AGB objects consistent with the influence of a binary? (6) Are fast out- 
flows contaminated by material that could represent accretion disk residue of 
shredded low mass companions? (7) Arc time scales of observed outflow pre- 
cession consistent with the gravitational influence of a binary on a disk? (8) 
Can double peaked line profiles be detected to identify accretion disks within 
the launch region? (9) Can shrouded novae outbursts from a M2 > Merit com- 
panion feeding the primary be detected in X-rays? (10) Improved statistics 
on the fraction of bipolar pPNe, the fraction of suitable precursor binaries for 
CEE, and the fraction of stars which evolve to be pPNe will constrain theories 
and assess whether all PNe incur asymmetry. 
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